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We consider the case of a density-driven metal-superconductor transition in the proximity of an 
electronic phase separation. In particular we investigate the interplay between superconducting 
fluctuations and density fluctuations, which become quantum critical when the electronic phase 
separation vanishes at zero temperature into a quantum critical point. In this situation the critical 
dynamical density fluctuations strongly affect the dynamics of the Cooper pair fluctuations, which 
acquire a more singular character with a z = 3 dynamical critical index. This gives rise to a scenario 
that possibly rules the disappearance of superconductivity when the electron density is reduced by 
elecrostatic gating at the LaAlOs/SrTiOs interface. 


PACS numbers: 74.40.Kb, 74.81.-g, 75.40.Gb, 74.62.-c 

I. INTRODUCTION 

There are many electronic systems where some kind 
of order is established upon changing the electron den¬ 
sity. Traditional examples are given by disordered doped 
semiconductors or MOSFETs, where a metal-insulator 
transition is found upon reducing the electron densit}^ 
and by heavy fermion systems, where competing phases 
(like, e.g., antiferromagnetism) are often spoiled by dop¬ 
ing, giving rise to metallic states often accompanied by 
superconductivit}El. High temperature superconducting 
cuprates and pnictides are other noticeable examples 
of systems where several electronic phases are realized 
by changing doping. More recently it was also found 
that the two-dimensional electron gas (2DEG) formed 
at the interface of transition-metal oxides (like, e.g., 
LaXOa/SrTiOs, with X=Al,Ti, LXO/STO) become su¬ 
perconducting when the electron density is varied by elec¬ 
trostatic gatin^^. A zero-temperature Quantum Griti- 
cal Point (SG-QGP) separates the superconducting phase 
from the non-superconducting on^. It may happen that 
the electronic liquid is not homogeneously distributed in 
the system, as in the presence of an electronic phase 
separation (PS) for instance. In this case, if the elec¬ 
tronic inhomogeneities are static and sizable, i.e., large 
enough to sustain the ordered phase (e.g., superconduc¬ 
tivity) one may expect that an inhomogeneous ordered 
state can be realized. Eor instance, one can envisage a 
state where superconducting ’’puddles” are embedded in 
a non-superconducting metallic phase of lower electronic 
density, and link together to establish an overall coher¬ 
ent state, which is destroyed when the electron density is 
r educed. 1^. 

The phase separation regime itself may be tuned by an 
external parameter, such as an electrostatic gate voltage, 
which controls the average electronic density, and dis¬ 
appear into a zero-temperature critical point (PS-QGP). 
Around this point the density fluctuations have a quan¬ 
tum nature and keep a dynamical character, giving rise 
to slow critical density fluctuations at very long wave¬ 


lengths. Of course, this physical situation can only oc¬ 
cur (at least approximately on large, even though not 
infinite, scales) if the Goulombic repulsion between elec¬ 
tron density fluctuations is substantially screened. In 
this case the critical electron density fluctuations in the 
clean limit dispW a typical behavior characterized by 

their propagatoi^El 

D(q, cj) = ((5n(q, Co’)(5n(—q, —cj)) 


where c„ is a fermionic scale, q= |q|, and m° = Cn^~'^ 
being the density correlation length) rules the dis¬ 
tance from the PS-QGP, where it vanishes. To be more 
specific, in LXO/STO this distance can be tuned by the 
gating potential Vg, oc {Vg —Vg), where Vg is the 
critical gating at which the QGP occurs. Near this QGP 
the Landau damping induces a term i\uj\/q. The largest 
contribution of these fluctuations (i.e. the maximum of 
their spectral density) occurs when the damping term is 
of the same order of the term. This leads to an over¬ 
damped dynamics with uj ^ q^ and therefore a dynamical 
critical index z =‘S characterizes the critical density fluc¬ 
tuations. 

Here, we address one single question: what happens if 
the QGP for phase separation occurs near the QGP of the 
density-driven ordering transition? Eor the sake of con¬ 
creteness, we will explicitly consider the case of density- 
driven superconductor-to-metal transition, but several of 
our arguments will apply to other transitions as well. 
The rather generic finding is that the most critical (i.e., 
singular) fluctuations imprint the dynamics of the less 
critical fluctuations. In the specific case of PS-QGP and 
SG-QGP proximity we will see that the z = 2 supercon¬ 
ducting (Gooper) fluctuations inherit the z = ‘5 dynamics 
of the underlying density fluctuations. 

In order to provide a first intuitive glimpse of the situa¬ 
tion, we first consider a toy model based on superconduc¬ 
tivity instantaneously taking place as soon as a (dynam¬ 
ical but slow) density fluctuation increases locally the 



2 


electron density above a critical value (Sect. II). In Sect. 
Ill we will provide general arguments to support the idea 
that the two QCP (the one closing the PS dome and the 
one related to density-driven superconductivity) tend to 
attract and coalesce. Sect. IV contains our concluding 
remarks. 


II. THE FLUCTUATING PUDDLE MODEL 

We start from a simple “toy model” model, where the 
critical character of the density-driven ordering is ne¬ 
glected, in order to make the above “imprinting” mech¬ 
anism more explicit. We assume that when a certain 
critical density ric is locally reached, a superconducting 
(SC) order parameter instantaneously arises. In partic¬ 
ular, it may happen that a dynamical density fluctua¬ 
tion is created, where locally the density is sufficiently 
high to sustain superconductivity. In this case, neglect¬ 
ing the transient behavior, one faces a superconducting 
puddle, which, however, only survives as long as the den¬ 
sity fluctuation survives. It is then natural to infer that 
the dynamics of the density fluctuation imprints (actu¬ 
ally it determines) the dynamics of the superconducting 
puddles and leads to an apparent ^ = 3 behavior for SC 
dynamics. 

Owing to the propagator in Eq. 0 - the density fluctu¬ 
ations in momentum space [i.e., the Fourier transform of 
Jn(r, t) = n(r, t) —n, n being the average density] satisfy 
the differential equation 

^dtSn{q, t) = + q^) (5n(q, t), ( 2 ) 

where = m^jcn^ and we take the momenta in units 
of Fermi wavevectors and time in units of inverse Fermi 
energy. Eq. 0 can be solved giving, for t > 0 , 



Figure 1. Time dependence of the puddle size in the case of 
superconducting puddles at any positive density {mPsc — 0)- 
The long dashed (red) straight line identifies the slope 1/3 
characteristic of the z — diffusive dynamics. The short 
dashed (green) straight line identifies the slope 1 characteris¬ 
tic of the z — 1 diffusive dynamics. The thick black curve is 
for = 10“^, while the thin blue curve is for = 10“^, 
in inverse of unit lengths square. Inset: spacial profile of the 
order parameter at various times, from t=10tot=10x2^^ 
(in units). The initial size of the fluctuating puddle is of 
3 unit lengths (in units). 


if the PS-QCP occurs at slightly lower density 
one needs a sufficiently large density fluctuation Sn > 
= m^c/d order to create the SC order in the 
region of the fluctuation. Thus, inside the region with 
(sufficiently) positive (5n(r, t) > g, one has 


[A(r,i)]^ = 


gSn{r,t) 


e {gSn{r, t) - m°^), (5) 


6n{q,t) = 5n{q,t = 0)e +« )*, ( 3 ) 

The corresponding real-space density can be coupled to 
the SC order parameter A(r)| by a term of the form 

HnA =-g J Sn{r)\A{r)fd‘^r, (4) 

while the SC order parameter itself is ruled by a Hamil¬ 
tonian 


Hsc = 


/ 


m°,|A(r)|2 + -|A(r)|4 


d\, 


where rules the distance from the SC critical point. 
Here space and time dynamics of the SC order parameter 
are purposely neglected, within our toy model, in the 
absence of coupling to the density fluctuations. If the 
SC-QCP and the PS-QCP coincide at the same critical 
density ric, and the SC phase is located at higher density, 
a positive fluctuation of Sn = n — ric > 0 instantaneously 
induces a static SC order parameter. On the other hand. 


where 0{x) is the Heaviside function. To test the occur¬ 
rence of a ^ = 3 dynamics for the SC order parameter, for 
the sake of definiteness, we consider the case when super¬ 
conductivity takes place in the puddle as soon as Jn > 0 
(i.e. = 0). We numerically generated at t = 0 in one 

dimension a density fluctuation with a Gaussian space 
profile and, after Fourier transforming, we let it evolve 
in time according to Eq. §. Transforming back to real 
space, we insert the new density profile in Eq. 0 and we 
identify the time evolution of the frontier of the SC pud¬ 
dle, obtaining the curves of Fig. One can see that, 
after a short time transient, a 2 : = 3 diffusive dynamics 
(the long-dashed line is a ^ = 1.85curve for compar¬ 
ison) occurs with the front of the SC puddle increasing 
like va ^ (black and blue curves corresponding to 
= 10“^ and = 10“^ respectively). The inset 
reports the profile of the SC order parameters A{x) at 
various times in the case = 10“^. In this case one 
can also see that, due to the finite extension of the crit¬ 
ical density fluctuation (^^ = 10 ^/^ unit lengths), a long 
time crossover to a ^ = 1 propagating behavior occurs 
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for t > 10^ (short-dashed straight line). 

A word is now in order to comment on the space and 
time units for systems like LXO/STO. The small densi¬ 
ties involved (n 0.01 — 0.08 el/cell) are such that the 
Fermi energies are small (a few tens of meV) and the 
inverse Fermi wavevector is roughly of order of 5-20 lat¬ 
tice spacings. Thus the extension of the puddle before 
the z = 1 crossover is reached can be estimated to be 
300-3000 lattice units for = 10“^ or = 10“^, re¬ 
spectively. For instance, in LTO/STO 300 lattice units 
correspond to about 120 nm. 


III. PERTURBATIVE HINTS 

In the above model one finds that a static SC order pa¬ 
rameter acquires a critical dynamics “riding” on a z = 3 
density fluctuations. We now analyze from a simple per¬ 
turbative point of view the interplay between the PS and 
the SC quantum criticality recovering the critical char¬ 
acter of the superconducting fluctuations, which was ne¬ 
glected in the above “toy model”. We consider the case 
when the superconducting order parameter fluctuates 
both in modulus and phase (like in the original Ginzburg- 
Landau case) and couples to nearly critical density fluctu¬ 
ations of an otherwise homogeneous metal. This scenario 
involves the interplay between standard Cooper fluctua¬ 
tions and critical density fluctuations. An obvious pre¬ 
requisite for this interplay is that the PS-QCP and the 
SC-QCP are in close proximity. This naturally raises the 
question whether such proximity requires a fine tuning of 
the model parameters or it naturally occurs due to some 
intrinsic “attraction” of the two critical point. The first 
result of our perturbative analysis is that the second pos¬ 
sibility is quite likely. First of all, it is well known that SC 
pairing is favored (and it may even be induced) by critical 
density fluctuations. It has been shown in the framework 
of strongly correlated electrons that near a PS-QCP su¬ 
perconducting pairing is induced due to th e mediation of 
the nearly critical density fluctnation j^^ l ^^ l On the other 
hand, it has also been proposecP^ that a density-driven 
superconducting transition may favor phase separation. 
Here we perturbatively show that the density fluctuations 
are made softer (and the compressibility correspondingly 
increases) when the SC-QCP is approached. We calculate 
the diagram (a) in Fig. representing the first pertur¬ 
bative corrections to in the D propagator of Eq.Q 
(which corresponds to the bare inverse compressibilit;^. 
The dashed lines represent the usual Cooper fluctuation 
propagatoi^ 

L(q,cj) = (A(q,cj)A(-q, -u)) 


cW-i\^\^rn^sc ^ ^ 

which clearly displays a z = 2 (cj ~ g^) overdamped 
dynamics. We find that, starting from separated QCPs, 
the finite bare inverse compressibility tends to be 


wvi; 


'VWh 




Figure 2. (Above) Self-energy correction of the density- 
density propagator (wavy line) due to the Cooper fluctuations 
(dashed lines), (below) Self-energy corrections of the Cooper 
propagator due to density fluctuations. 


reduced by the coupling to the superconducting critical 
fluctuations 


0 2 1 ^ 
mn = m^- g log —^ 

yielding a smaller dressed inverse compressibility rUn- 
Here and also measure the distances from the 
PS-QCP and SC-QCP respectively, while A is a non- 
universal high-energy cutoff. Therefore, no matter how 
weak is the coupling between superconducting and den¬ 
sity fluctuations, the inverse compressibility = rrin 
tends to decrease when approaching the SC-QCP (where 
m^c ^0)- A pictorial view of this scenario is reported 
in Fig. Although a full renormalization group calcu¬ 
lation is in order to draw definite conclusions, the above 
perturbative finding indicates that the proximity of the 
PS-QCP and the SC-QCP is not a fortuitous (and un¬ 
likely) coincidence. 

Based on the above indication, we now investigate the 
interplay of the SC and density fluctuations assuming 
that the the SC and PS criticality are sufficiently close 
to strongly influence each other and we impose a formal 
coincidence of the two QCPs. Thus we carry out pertur¬ 
bative calculations in two dimensions starting from the 
standard Cooper fluctuation propagator in Eq. (|^, with 
z = 2, and the density fluctuation propagator of Eq. Q, 
with z = 3. To impose the QCP coincidence maximiz¬ 
ing the effect of the critical interplay we set to zero both 
and in the calculation of the diagram in Eig. 
(b). We find that the SC fluctuations are strongly renor¬ 
malized and at first order their propagator acquires the 
form 


Z(q,a;) 


1 

CscQ - ico^/q^ ' 


Again the largest contribution of these fluctuations (i.e. 
the maximum of their spectral density) occurs when 
the damping term is of the same order of the q 

term leading to an overdamped dynamics with ^ q^. 
And again this corresponds to a dynamical critical index 
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Figure 3. (a) Schematic phase diagram of the density-driven 
superconductivity with a SC-QCP at a critical density ric- 
(b) Chemical potential vs. density displaying an increase of 
the compressibility Hi (dashed line) in correspondence to the 
SC-QCP 


2 : = 3. Notice that the constant terms in the denomi¬ 
nator, which would drive the system away from the crit¬ 
ical point has been self-consistently set to zero to keep 
the criticality condition fulfilled. Once more we stress 
that a full renormalization group treatment is required 
to firmly establish the behavior near the two QCP, but 
our perturbative calculations are clearly indicative of the 
following scenario. The two QCP tend to coalesce and 
would generically give rise to a first-order transition with 
a PS between a low-density metallic state and a high- 
density SC state. However, as it usually occurs, a tuning 
of some parameter (like gating or magnetic field) may 
bring the system near a quantum tricritical point, where 
the SC critical temperature vanishes and, at the same 
time, the compressibility diverges, n = m~^ oc (i.e., 
the two QCPs merge). At this point the SC fluctuations 
acquire a 2 ; = 3 dynamics, while the effective dimension¬ 
ality d + 2 ; = 5 is larger than the upper critical dimension 
due which is three for a tricritical point. In this case, for 
a two-dimensional system, one would expect mean-field 
critical indexes. 


Figure 4. Density n vs. gating Vg phase diagram. The SC 
critical temperature Tc(n) vanishes along the line of QCP 
represented by the black dotted line. In the PS region static 
puddles of denser superconducting phase are embedded in a 
less dense normal metallic phase. The darker shaded area rep¬ 
resents the region where PS-QCP interplay with the nearby 
SC-QCPs giving rise to mixed dynamical quantum SC and 
density fluctuations. 


IV. POSSIBLE REALIZATION: LXO/STO 
INTERFACE 

All the above analysis is based on the presence of 
long-ranged nearly critical density fluctuations. Obvi¬ 
ously such fluctuations are only possible if the long- 
range Coulombic (LRC) electron-electron interaction is 
strongly screened, otherwise these fluctuations would be¬ 
come high-energy barely damped plasmonic exci tation s. 
In this case the dynamical critical index 2 ; = l.CSISl R 
has also been showrP^ that a moderate screening, like 
the one due to a metallic gate close to a Josephson- 
Junction array, may still lead to 2 : = 1 and a rather large 
value of the correlation length critical index v > 

The question therefore becomes: to what extent our 
strong-screening assumption is valid leading to the strong 
Landau-damping effects responsible for the large, 2 : = 3, 
value of the dynamical critical index? The answer nat¬ 
urally depends on the physical system at hand. If the 
SC transition occurs in a very low-density system, LRC 
interaction likely plays a relevant role and the scenario of 
the previous sections would hardly apply. On the other 
hand, when the SC transitions takes place in a metallic 
state with high mobility and relatively high electron den¬ 
sities, all interactions become short ranged and the super¬ 
conducting fluctuations may be strongly damped. The 
text-book case of the metal-to-superconductor transition 
driven by Cooper fluctuations naturally realizes this case 
with z = 2 [cf. Eq. §]. The same critical index charac¬ 
terizes the superconducting QCP of Ref. [191 LXO/STO 
interfaces, where n ~ 10^^ — 10^^ electrons/cm^ likely 
belong to this class systems where the screening plays a 
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relevant role thereby leaving the possibility of large criti¬ 
cal density fluctuations open. Moreover, one should also 
remember that SrTiOs, where the interface charges re¬ 
side, has a huge dielectric constant ranging from 300 to 
25000 strongly enhancing the screening. 

Disorder is another relevant issue to be discussed. We 
consider systems where the microscopic disorder is weak 
and rather homogeneously distributed. In this case the 
metallic character of the systems is poorly affected by the 
impurity scattering at most resulting in weak localization 
of no relevance to our purposes (like in LTO/STO). On 
the other hand, the presence of a PS naturally entails 
the occurrence of a strongly inhomogeneous state emerg¬ 
ing from this otherwise weakly disordered metal. This 
inhomogeneity would poorly affect transport properties 
because it mostly enhances forward scattering leaving the 
mobility high. However such static inhomogeneity would 
greatly affect the critical properties raising the question 
of the Harris criterion to be satisfied in order to have a 
true phase transition. The metallic state near the PS- 
QCP is not disordered because the density fluctuations 
are slow (that is critical), but fully dynamical excitations 
of the electron gas and they are not a static (quenched) 
source of disorder. In this case the Harris criterion simply 
does not apply (at least, in its standard form). In partic¬ 
ular, when the metal-SC transition is driven by the gating 
potential a = 3/2 scaling is observed.l^In this case 


one can argue that gating brings the system close to a 
PS-QCP, where no static puddles are present (thus Har¬ 
ris criterion is not applicable) and the density critically 
fluctuates in a clean metallic state. In this scheme the 
experimental observation of = 3/2 could arise accord¬ 
ing to the scenario of Sect. |m] from a z = 3 dynamical 
critical behavior of the superconducting fluctuations to¬ 
gether with a mean-field like exponent v = 1/2. We find 
it appealing to consider and explore the possibility of this 
alternative scenario, in which the superconducting fluc¬ 
tuations are coupled to (and imprinted by) nearly critical 
density fluctuations. 

We finally remark that our scenario may also apply 
to those cases where superconducting fluctuations cou¬ 
ple to longitudinal nematic modes (also characterized by 
z = possibly explaining the occurrence of a super¬ 
conducting QCP with zu = 3/2 reported in a gate doped 
La 2 -ccSra,Cu 04 sample.!^ 
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